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Glossary 
Cellular Automata 

Discrete variables on a discrete lattice change in discrete time steps. 
Ising model 

Neighbouring variables prefer to be the same but exceptions are possible. 
The probability for such exceptions is an exponential function of "tempera- 
ture" . 

Percolation 

Each site is randomly either occupied or empty, leading to random clus- 
ters. At the percolation threshold for the first time an infinite cluster is 
formed. 

Universality 

Certain properties are the same for a whole set of models or of real objects. 
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1 Definition 



This article introduces into the whole section on Social Sciences, edited by 
A. Nowak for this Encyclopedia, concentrating on the applications of mathe- 
matics and physics. Here under "mathematics" we include also all computer 
simulations if they are not taken from physics, while physics applications 
include simulations of models which basically existed already in physics be- 
fore they were applied to social simulations. Thus obviously there is no 
sharp border between applications from physics and from mathematics in 
the sense of our definition. Also social science is not defined precisely. We 
will include some economics as well as some linguistics, but not social insects 
or fish swarms, nor human epidemics or demography. Also, we mention not 
only this section by also the section on agent-based modelling edited by F. 
Castiglione as containing articles of social interest. 

2 Introduction 

If mathematical/physical methods are applied to social sciences, a major 
problem is the mutual lack of literature knowledge. Take for example the 
Schelling model of racial segregation in cities Sociologist don't cite the 
better and simpler Ising model, physicists ignored the Schelling model for 
decades, and sociologists also ignored better sociology work [2]. For sim- 
ulations of financial markets, many econophysicists thought that they had 
introduced Monte Carlo and agent-based simulations to finance, not knowing 
of earlier work from some forward-looking Nobel laureates in economics |3l H] . 
For inter-community relations, already 25 centuries ago analogies with liquids 
were pointed out by Empedokles in Sicily [S]. More recently, Ettore Majo- 
rana [6] around 1940 suggested to apply quantum-mechanical uncertainty 
to socio-economic questions. With emphasis shifted to statistical physics, 
sociophysics and econophysics became fashionable around the change of the 
millennium, but continuous lines of research by some physicists started [7] 
already in 1971. In the same year Journal of Mathematical Sociology started 
and published Schelling's model of urban segregation [1], which is a modifi- 
cation of the Ising magnet at zero temperature. 1982 saw the start of two 
other lines of research by physicists on socio-economic questions [8[[9]. 

Languages have been simulated on computers by decades, while the in- 
terest of physicists is more recent [101 [H], triggered mostly by a model of 
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language competition fn\ . 

We do not mention chemists since at present they play no major role 
in this field. However, the 1921 chemistry Nobel laureate F. Soddy [13], to 
whom we owe the "isotope" concept, already worked on economic, social and 
political theories, and his finance work of the 1930's was still cited in 2007. 
The present authors try it the other way around: First apply physics to social 
sciences, and then get the Nobel prize (for literature: science fiction.) 

3 Some Models and Concepts 

Physicist Albert Einstein said that models should be as simple as possible, 
but not simpler. In this spirit we now introduce some basic physics models 
and concepts for readers from social sciences. They don't have to study 
physics for many years, the following examples give the spirit. All models are 
complex in the sense that the behaviour of large systems cannot be predicted 
from the properties of the single element. 

3.1 Cellular Automata 

Mathematicians denote cellular automata often as "interacting particle sys- 
tems", but since many other models or methods in physics use interacting 
particles, we do not use this term here. A large c?-dimensional lattice of L'^ 
sites carries variables Si {i = 1,2, . . . L'^) which can be either zero or one; 
more generally, they are small integers between 1 and Q > 2. The lattice 
may be square (four nearest neighbours), triangular (six nearest neighbours), 
or simple cubic (also six nearest neighbours, but in c? = 3 dimensions); many 
other choices are also possible. Time t = 1,2, ... increases in steps. At each 
time step, each ^^(t + l) is calculated anew from a deterministic or probabilis- 
tic rule involving the neighbouring Sk{t) of the previous time step. This way 
of updating is called "simultaneous" or "parallel"; one may also use sequen- 
tial updating where Si depends in the current values of Sk] then the order of 
updating is important: random sequential, or regular like a typewriter. 

An example is a biological infection process: Each site i becomes perma- 
nently infected, S'j = 1, if at least one of its nearest neighbours is already 
infected. Computers handle that efficiently if each computer word of, say, 32 
bits stores 32 sites, and if then 32 possible infections are treated at once by 
bit-by-bit logical-OR operations [H]. 
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3.2 Temperature 



We know temperature T from the weather reports, but in physics it enters 
according to Boltzmann into the probabihty 



to observe some configuration with an energy E. Here T is the temperature 
measured in Kelvin (about 273 + the Celsius or centigrade temperature), and 
ks the Boltzmann constant relating the scales of energy and temperature. 
For simplicity we now set — 1, i.e. we measure temperature and energy 
in the same units. If g different configurations have the same energy, then 
S = ln{g) is called the entropy, and the probability to observe this energy is 
oc g exp{—E /T) = exp(— F/T) with the "free energy" F = E — TS. 

In a social application we may think of peer pressure or herding: If your 
neighbours drink Pepsi Cola, they influence you to also drink Pepsi, even 
though at present you drink Coca Cola. Thus let E be the number of nearest 
neighbours drinking Pepsi Cola, minus the number of Coke drinking nearest 
neighbours. The probability for you to switch then is given by the energy 
difference and equal to exp(— 2£^/T) (or 1 if £^ < 0) in the Metropolis algo- 
rithm, or 1/(1 + exp(2ii^/T) in the Glauber or Heat Bath algorithm. In both 
cases there is a tendency to decrease E. In the limit T — one never makes 
a change which increases E, while for small positive T one increases E with 
a low but finite probability. In the opposite limit of infinite temperature, the 
energy becomes unimportant and all possible configurations become equally 
probable. Neither zero nor infinite temperature are usually realistic. 

In this sense, decreasing the energy E is the most simple or most plausible 
choice, and the temperature measures the willingness or ability to deviate 
from this simplest option, e.g. to withstand peer pressure. But temperature 
also incorporates all those random accidents of life which infiuence us but 
are not part of the social model. For example, it may happen that there is 
no Pepsi Cola available even though all your neighbours drink Pepsi and you 
want to follow them. Investors have to make their financial choices under the 
infiuence of their clients, whose life is shaped by births, marriages, deaths, 
or other personal events which are not included explicitly into a financial 
market model. These accidents are then simulated by a finite temperature, 
entering the probability that one does not follow the usual rule. 

The ability to withstand peer pressure and the randomness of personal 
lives are in principle two different things, and if one wants to include them 
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both one needs two different temperatures Ti and T2, which do not exist in 
traditional physics [T^ . 

3.3 Ising Model 

In the model published by Ernst Ising in 1925, the variables Si are not or 
1, but ±1: 

E = -J2S,Sk-Bj2S, (2) 

i,k i 

and for i? = this corresponds to the above Coke versus Pepsi example. 
The first summation runs over all neighbour pairs, the second over all sites. 
Thus if site i considers changing its variable, the energy change is ±AE = 
2{J2k — B) and enters through exp(— Aii^/T) into the probabilities to flip 
Si] now k runs over the neighbours of i only. (If instead of flipping one Si 
one wants to exchange two different variables Si and Sj, moving Si into site 
i and Sj into site then one has to calculate the energy changes for both 
sites i and j in this "Kawasaki" kinetics.) A computer program and pictures 
from its application are given elsewhere in this Encyclopedia 

In physics, the Si are magnetic dipole moments of the atoms, often called 
spins, and B is proportional to the magnetic field. Usually, physicists write an 
exchange constant J before the first sum, but we set J = 1 for simplicity here. 
The model was invented to describe ferromagnetism, like in the elements iron, 
cobalt or nickel. Later it was found to describe liquid-vapour equilibria and 
other phase transitions. We know that iron at room temperature is magnetic, 
and this corresponds to the fact that for < T < Tc and zero field B the 
Ising model has the majority of its spins in one direction (either mostly +1 or 
mostly -1), while for T > half of the spins point in one and the other half 
in the opposite direction. The magnetisation M = J^i^i, often normalised 
by the number L'^ of spins, is therefore an order parameter. The critical 
temperature Tc is often named after Pierre Curie. 

In one dimension, we have Tc = 0; in the square lattice in two dimensions 
we know Tc = 2/ ln(l + ^/2) exactly, while on the simple cubic lattice Tc ~ 
4.5115 is estimated only numerically. Of course, one has generalized the 
model to more than nearest neighbours, to more than two states ±1 for each 
spin, and to disordered lattices and networks. 
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3.4 Percolation 



Simpler than the Ising model but less useful is percolation theory, reviewed 
more thoroughly in this Encyclopedia in the section edited by M. Sahimi. 
Each site of a large lattice is randomly occupied with probability p, empty 
with probability 1 — p, and clusters arc sets of occupied neighbouring sites. 
There is one percolation threshold Pc such that for p < Pc only finite clusters 
exist, for p > Pc also one infinite cluster, and at p = Pc even several infi- 
nite clusters may co-exist, which are fractal: The number of occupied sites 
belonging to the infinite clusters varies at Pc as where D is the fractal 
dimension. Here "infinite" means: spanning from one end of the sample of 
L'^ sites to the opposite end, or: increasing in average number of sites with 
a positive power of L. In one dimension, again one has no phase transition 
(pc = 1), on the square lattice Pc ^ 0.5927462 and on the simple cubic lattice 
Pc ^ 0.311608 are known only numerically, with a fractal dimension of 1, 
91/48 and ^ 2.5 in one to three dimensions. 

In the resulting disordered lattices, each site has from to z neighbours, 
where z is the number of neighbours in the ordered lattice p — 1. If one 
neglects the possibility of cyclic links one finds Pc — l/{z — 1) in this Bethe 
lattice or Cayley tree. Near this percolation threshold the critical exponents 
with which several quantities diverge or vanish are the same as in the random 
graphs of Erdos and Renyi. But this percolation theory was published nearly 
two decades earher, in 1941, by the later chemistry Nobel laureate P. Flory. 

3.5 Mean Field Approximations 

What is called "mean field" is called "representative agent" theory in eco- 
nomics, and is widespread in chemistry where the changes in the concentra- 
tions of various reacting compounds are approximated as functions of these 
time-dependent concentrations. A particularly simple example is Verhulst's 
logistic equation dx/dt = ax{l — x) , known as Bass diffusion in economics. 
We now explain why this approximation is unreliable. 

Let us return to the above Ising model of Eq.(2) and replace the Sk there 
by it's average < Sk >— m — M/L'^ which is a real number between -1 and 
+1 instead of being just -1 or -|-1; m is the normalised magnetisation. Then 
the total energy E is approximated as the sum over single energies Ei. 

E = Y,Ei Ei^{-J2<Sk> -B)Si = -B'Si 

i k 
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with a mean magnetic field B' = B + J2k < >= B + mz where z again 
is the number of lattice neighbours. The system now behaves as if each spin 
Si is in an effective field B' influenced only by the average magnetisation m 
and no longer directly by its neighbours 5*^. The two possible orientations of 
Si have the energies ±5', giving an average 

m =< Si >= tanh{B'/T) = tanh[(5 + zm) /T] (3a) 

and thus a self-consistency equation for m. Expanding the hyperbolic tangent 
into a Taylor series for small m and B we get 

5 = (1 - z/T)m + + ... (36) 

which gives a Curie temperature = z, since for T < Tc the magnetisation is 
m = ±[3(z/T — 1)]^/^ oc (Tc — T)^/^. Similar approximations for liquid- vapour 
equilibria lead to the Van der Waals equation of 1872, which may be regarded 
as the first quantitative theory of a complex phenomenon, (m there is the 
difference between the liquid and the vapour density.) Nowhere in Eqs(2,3a) 
have we put in that there is a phase transition to ferromagnetism; it just 
arises from the very simple interaction energy SiSk between neighbouring 
spins, and similarly the formation of raindrops emerges from the interaction 
between the molecules of water vapour. The water molecule is the same H2O 
in the vapour, the liquid or the ice phase. 

But this nice approximation contradicts the results mentioned above. For 
the chain, square and simple cubic lattice it predicts a T^, = ^ = 2, 4 and 6 
while the correct values are 0, 2.2, and 4.5. Particularly in one dimension 
it predicts a phase transition at a positive Tc while no such transition is 
possible: Tc = 0. This was the main result of Ernst Ising's thesis in 1925. 
And even in three dimensions, where the difference in Tc between 4.5 and 6 
is less drastic, the above square-root law for m is wrong, since m varies for 
T slightly below Tc roughly as (Tc — T)°''^^. Thus mean field theory. Van der 
Waals equation, and similar approximations averaging over many particles 
are at best qualitatively correct. They become exact when each particle 
interacts equally with all other particles. 

Analogously for percolation, Flory's approximation of neglecting cyclic 
links and the Erdos-Renyi random graphs lead to results corresponding to 
mean field approximations and should not be relied upon in one, two, or 
three dimensions with links between nearest neighbours only. 

For cellular automata a particularly drastic failure of analogous mean 
field approximations (differential equations) was given by Shnerb et al [TB] 
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for a biological problem. Even simpler, many cellular automata on the square 
lattice lead to blinking pairs of next-nearest neighbours: at even times one 
site of the pair is 1 and the other is 0, while at odd times the first is and 
the second is 1. Averaging over many sites destroys these local correlations 
which keep the blinking pair alive. 

4 Applications 

A thorough review of "sociophysics" was given recently by Castellano et al 
[T7] . a long list of reference by Carbone et al ^8]. Some work of social 
scientists is reviewed by Davidsson and Verhagen in the section on agent- 
based simulations in sociology, while Troitzsch in this section reviews both 
social scientists and physicists. His book with Gilbert [29j is, of course, more 
complete. Thus we merely sketch here some the areas covered in greater 
detail in the other articles or in the cited literature. 

4.1 Elections 

A social scientist may be interested to predict the fate of one particular party 
or candidate in one particular election, or to explain it after this election. 
A physicist, accustomed to electrons, hydrogen atoms and water molecules 
being the same all over the world may be more interested to find which 
universal properties all elections have in common. Figure 1, kindly sent 
by Santo Fortunato, is an example. Let v be the number of votes which 
a candidate got, Q the number of candidates in that election, and the 
total number of votes cast. Then the probability distribution P{v, Q, N) 
for the number of votes is actually a function f{vQ/N) of only one scaled 
variable, and that variable vQ/N is the ratio of the actual number v of 
votes to the average number N/Q of votes per candidate. Various countries 
and various electsions, all using a proportional election system, gave the 
same curve f{vQ/N) which is a parabola on this double-logarithmic plot and 
thus corresponds to a log-normal distribution. In Brazil, however, where the 
personality of a candidate play a major role, not only the party membership, 
the results were different. These authors [19] also present a model to explain 
the log-normal distribution. 

Other models for opinion dynamics are reviewed elsewhere in this Ency- 
clopedia 0. 
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Santo Fortunate & Claudio Castellano, PRL 99, 138701 




vQ/N 



Figure 1: The vote distribution in several countries and elections is a function 
only of the scaled variable vQ/N . From |19j . 

4.2 Financial Markets 

Agent-based simulation of stock markets [2D] are a typical example of com- 
plex systems applications: In these models not the single agent but their 
(unconscious) cooperation produces the ups and downs on the stock market, 
the bubbles and the crashes. These models deal with the more or less random 
fluctuations, not with well founded market changes due to new inventions or 
major natural catastrophes. 

Real markets give at each time interval a return r which is the relative 
change of the price. Typically, an index of the whole market like Dow Jones 
changes each trading day by about one percent. Much larger fluctuations 
are more rare, and the probability to have a change larger than r decays 
for large r as Fat tails. The sign of the change is barely predictable, 

but its absolute value is: Volatility clustering. Thus in calm times when |r| 
was small, tomorrow's |r| probably is also small, whereas for turbulent times 
with high |r| in the past one should also expect a large |r| tomorrow. The 
daily weather behaves similarly: presumably tomorrow will be like today. 
Perhaps even multifractality exists in real markets, similar to hydrodynamic 
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201 * 201 lattice witti concentrations summed from 1 to 59 percent; and Zhang square root relation 




chiange 



Figure 2: Simulated return distribution in the Cont-Bouchaud percolation 
model of stock markets ^21]. The asymptotic slope to the right is about —2.9. 



turbulence. 

A simple model, going back to Bachelier more than a century ago, would 
throw a coin to determine whether the market tomorrow will go up or down. 
This simple random-walk or diffusion model was shown by Mandelbrot in the 
1960's not to describe a real market; it lacks fat tails and volatility clustering 
but may be good for monthly changes. Many better agent-based models have 
been invented during the last decade and reproduce these real properties. Fig. 
2; the Cont-Bouchaud model is based on the above percolation theory [21] . 
while the Minority Game tells you it is better not to be with the big crowd 
[221. 



4.3 Languages 

The versatility of human languages distinguishes us from the simpler commu- 
nication systems of other living beings. With computers or mathematically 
exact solutions [23] models have been studied for the learning of a language 
by children or for the evolution of human languages out of simpler forms. 

Closer to simulations in biology with the Darwinian selection of the fittest 
are the models of competition between various languages of adult humans: 
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number of speakers 

Figure 3: Simulated size distributions for human languages (full circles, and 
line), compared with reality (open circles). From [25] . 

Will the Welsh language survive against English in Great Britain? Similar 
to Lotka-Volterra equations for prey and predator in biology, some nonlinear 
differential equations [12] seem to describe the extinction of the weaker lan- 
guage. Better statistics are available for the size distribution of languages, 
where " size" is the number of people speaking this language. Here one model 
of de Oliveira et al found good agreement with reality, Fig.3; other models 
[Mj were less successful, inspite of many simulations from physicists. 

5 Future Directions 

The future should see more work in what we have shown here through our 
three figures: Searching for universal properties, or the lack of them, in the 
multitudes of models and in reality. Biology became a real science when the 
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various living beings were classified into horses, mammals, vertebrates etc. 
Within each such taxonomic set all animals have certain things in common, 
which animals in other taxonomic sets do not share. This check for univer- 
sality is different from improving our ability to ride horses. Thus making 
money on the stock market, or explaining the crash of 1987, is nice, but 
investigating the exponents of the fat tails, Fig.2, of all markets may give us 
more insight into what drives a market and what differences exist between 
different markets. Winning one particular election and predicting the winner 
is important, but universal scaling properties as in Fig.l may help us to un- 
derstand democracy better. Preventing the extinction of French language in 
Canada is important for the people there, but explaining the overall statistics 
of languages in Fig. 3 is relevant globally. 

It is in these general aspects where the methods of mathematics and 
physics seem to be most fruitful. One specific problem is better solved by the 
local people who know that problem best, not by general simplified models. 
A useful future approach for interacting agents would be their realisation by 
neural network models [2B|. 
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